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SUMMARY 


A periodic-average  flow  measurement  technique  involving  a 
hot-wire  anemometer  system  was  used  to  measure  the  periodically 
unsteady  and  three-dimensional  fluid  velocity  field  between  blade  rows 
in  the  first  stage  of  a low-speed,  multistage,  axial-flow  research 
compressor.  These  data  suggest  that  the  fluid  flow  through  the 
imbedded  rotor  and  stator  rows  are  appreciably  unsteady,  in  a periodic 
fashion,  in  portions  of  the  compressor  annulus.  Illustrative  examples 
of  periodic-average  fluid  flow  field  variation  with  rotor  blade 
sampling  position  in  stop-action  sequence  are  presented  for  different 
locations  in  the  compressor.  A simple,  first  order  approximation 
physical  description  of  the  blade  wake  flow  transport  and  interaction 
process  largely  based  on  experimental  data  interpretation  is  proposed 
to  organize  and  to  help  explain  the  observations  made.  Blade  span 
variations  of  flow  data  reflect  end-wall  effects.  Inlet  guide  vane 
exit  flow  data  involve  some  unusual  unsteady  flow  effects. 
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radial  flow  angle  (Figure  10;  Equation  B-13),  degrees 

absolute  tangential  flow  angle  with  respect  to  axial 
direction  (Figure  10;  Equation  B-12) , degrees 
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1 . INTRODUCTION 


Modern  design  systems  for  turbomachines  rely  on,  to  a great 
extent,  empirical  correlations  to  reflect  real  fluid  flow  effects. 

This  state  of  design  competence  often  leads  to  good  machines. 

However,  more  often  than  can  be  afforded,  the  resultant  hardware 
involves  \mdesirable  deficiencies  that  can  be  related  to  Insufficient 
knowledge  about  the  fluid  mechanics  involved.  For  example,  mis- 
understanding of  the  substantial  relationships  between  the  unsteadiness 
of  turbomachine  flow  and  the  efficiency,  aerodynamic  and  aeroelastic 
stability,  and  noise  production  of  such  machines  can  result  in  costly 
disappointments.  Mikolajczak  [1]  suggests  that  further  significant 
improvements  in  turbomachine  technology  should  be  sought  through 
better  understanding  and  control  of  the  unsteady  flows  involved. 

Progress  in  fluid  flow  measurement  has  resulted  in  the  develop- 
ment of  a variety  of  techniques  (see,  for  example.  References  2-27) 
for  observing  the  unsteady  aspects  of  turbomachine  flows.  Precise 
coordination  of  data  acquisition  with  rotor  sampling  position  has  made 
possible  the  extraction  of  the  periodically  unsteady  flow  (periodic- 
average  flow)  data  from  the  entire  collection  of  information.  The 
periodic-average  flow  is  important  because  of  the  role  it  plays  with 
respect  to  forced  blade  vibration,  discrete  frequency  noise  generation, 
and  turbomachine  energy  transfer.  Data  yielding  detailed  information 
about  the  sequential  variation  of  periodic-average  flow  with  rotor 
sampling  position  are  rare  [14,15,16,24,26]. 

In  this  report  are  offered  detailed  periodic-average, three- 


I 


2 


I dimensional  flow  data  for  the  first  stage  of  a low-speed,  multistage, 

axial-flow  research  compressor  that  demonstrate  flow  field  variation 
with  sequential  change  in  rotor  sampling  position.  A physical 

I description  of  the  complicated  fluid  mechanical  processes  involved 

t 


is  proposed. 


i 


i 

t 2.  RESEARCH  COMPRESSOR  FACILITY 


The  research  compressor  facility  of  the  Iowa  State  University 
Engineering  Research  Institute/Mechanical  Engineering  Department 
Turbomachinery  Components  Research  Laboratory  was  used  for  the 
present  study-  Briefly  reviewed  below  is  the  research  compressor  and 
related  equipment  and  instrumentation.  For  more  detailed  information 
see  Reference  28. 


Axial  Flow  Research  Compressor 


Figure  1 is  a sketch  of  the  entire  research  conpressor  apparatus. 
Figure  2 depicts,  in  more  detail,  the  compressor  portion  only.  The 
inlet  guide  vane  (IGV)  row  and  three  identical  rotor-stator  stages 
were  within  an  annulus  having  constant  hub  (0.284  m)  and  tip  (0.406  m) 
diameters.  All  stationary  blade  rows  had  the  same  number  of  blades 
(37)  and  were  mounted  on  separate  ring  assemblies  which  could  be 
moved  independently  or  in  desired  combinations.  The  rotor  blade  rows 
each  included  the  same  number  of  blades  (38)  and  were  assembled 
together  so  that  all  of  the  respective  blade  stacking  axes  were  aligned 
axially.  All  of  the  blades  were  constructed  of  a plastic  material 
(Monsanto  ABS)  with  British  C4  sections  reflecting  a free  vortex  design. 
Blade  characteristics  are  summarized  below: 

Number  of  blades  per  row  IGV  and  stator  rows  - 37 

rotor  rows  - 38 

Blade  span  (constant)  6.10  cm  (2.4  in.) 


Blade  chord  (constant),  c 


3.05  cm  (1.2  in.) 
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Blade  section  maximum  thickness/  lOX 

chord  ratio,  t /c 
max 

Blade  geometry  details  are  tabulated  in  Table  1 with  associated 
variables  defined  in  Figure  3.  Measuring  stations  in  the  compressor 
were  positioned  "in  line"  axially  and  spaced  apart  as  shown  in  Figure  4. 

2.2.  Stationary  Blade-Row  and  Probe  Actuators 

Individual  or  combined  stationary  blade-row  motion  was  accomplish- 
ed with  the  circumferential  actuating  system.  The  actuator  consisted 
of  a rack  and  gearmotor-drlven  circular-arc  dovetailed  slider  moving 
within  a sturdy  base.  Slider  travel  was  monitored  with  linear 
potentiometer  output  voltage  correlated  to  motion  by  a linear  least 
squares  fit.  In  this  fashion,  the  slider  motion  could  be  adjusted  to 
within  0.05  . The  slider  was  connected  to  the  movable  IGV  and  stator 
blade  row  rings  with  adjustable  linkages  so  that  a variety  of  stationary 
blade  row  positioning  schedules  could  be  achieved.  Scales  on  each 
blade-row  ring  and  the  outer  compressor  casing  were  used  for  ascertaining 
the  precise  location  of  each  blade  row.  The  scales  were  calibrated  in 
terms  of  degrees,  positive  in  the  direction  of  rotation.  A reading  of 
0.0°  corresponded  to  the  stacking  axis  of  a predetermined  reference 
blade  in  each  row  being  "in  line"  with  the  measurement  stations  of  the 
compressor.  Any  circumferential  distance  of  a reference  blade  stacking 
axis  from  this  zero  location  was  denoted  as  YO,  its  circumferential 
position.  Throughout  these  tests,  IGV  and  stator  blade  row  rings 
were  positioned  for  minimum  sound  as  follows: 


9 


Inlet  Guide  Vanes 

''°IGV'=S 

= 0.00 

First  Stator 

= 0.17 

Second  Stator 

™2s'"s 

= 0.56 

Third  Stator 

™3s'=S 

= 0.77 

The  relative  positioning  of  some  of  these  blades  may  be  seen  in  Figure  5. 
As  will  be  explained  later,  the  circumferential  actuator  was  the  key 
to  providing  the  means  for  measuring  the  blade-to-blade  variations 
of  flow  in  the  compressor.  Since  the  probes  used  could  not  be 
actually  traversed  in  the  circumferential  direction,  the  stationary 
and  "sampled"  rotor  blades  had  to  be  moved  circumferentially  relative 
to  a fixed  probe. 

The  probe  actuator  (L.  C.  Smith  Con^jany  model  6180)  was  used 
with  the  control  indicator  (L.  C.  Smith  Company  model  DI-3R)  and 
actuator  switch  box  (L.  C.  Smith  Company  model  DI-4R-SB)  to  vary  the 
probe  yaw  angle  and  immersion  depth.  Probe  angles  and  positions  in 
the  compressor  annulus  were  monitored  by  observing  mechanical  digital 
counter  readings  and  linear  potentiometer  output  voltages.  Each 
potentiometer's  output  voltage  was  correlated  with  motion  with  a linear 
least  squares  fit  which  allowed  probe  angles  and  immersion  positions 
to  be  measured  within  0.05°  and  0.15  mm  respectively. 

2.3.  Pressure  and  Temperature  Measurement  Instrumentation 

All  working  fluid  pressure  measurements  were  made  with  precision 
water- in-glass  manometers  (Meriam  type  Incl. ) which  were  calibrated 
with  a micromanometer  (Meriam  type  Micro.).  Room  air  temperature  was 


Figure  5.  Blade  cascade  showing  relative  positions  of  blades  for 
several  rotor  sampling  positions. 
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measured  with  mercury-in-glass  thermometers,  while  working  fluid 
temperature  was  measured  with  a copper-cons tantan  thermocouple  and 
a precision  millivolt  potentiometer  (Leeds  and  Northrop  model  8686). 
Barometric  pressure  was  measured  using  a mercury-in-glass  barometer 
(Princo  Instruments,  Inc.  model  B-222) . 

2.4.  Per iodic- Average  Measurement  System 

The  periodic-average  measurement  system,  shown  schematically 
in  Figure  6,  was  composed  of  the  following  components: 

(1)  Single  slanted  hot-wire  probe  (Disa  model  55P02  Modified) 

(2)  Constant  temperature  anemometer  (Thermo-Systems,  Inc. 

[TSI]  model  lOlOA) 

(3)  Linearizer  (TSI  model  1072) 

(4)  Periodic  sample-and-hold  circuit 

(5)  Photoelectric  triggering  circuit 

(6)  Signal  averaging  circuit 

(7)  Digital  scanning  voltmeter  (Hewlett-Packard  model  3480D) 

(8)  Desk-top  calculator  (Hewlett-Packard  model  9821A) 

(9)  Oscilloscopes  (Tektronix  Inc.) 

The  hot-wire  probe  involved  a 5 pm  diameter  platinum-plated 
tungsten  wire  with  a 1.25  mm  sensing  portion  centered  on  the  probe  axis. 
The  wire  was  copper  and  gold  plated  at  the  ends  and  positioned  at 
an  angle  of  54.7°  from  the  probe  axis.  The  hot-wire  anemometer  and 
linearizer  were  used  to  produce  a linear  relationship  between  air 
velocity  and  output  voltage.  The  periodic  sample-and-hold  circuit  was 
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phase  locked  with  the  compressor  rotor  by  the  photoelectric  pickup 
which  was  triggered  by  a slot  in  a disc  rotating  with  the  compressor 
£ shaft.  The  photoelectric  pickup  was  connected  to  the  stationary 

blade  row  actuator  with  an  adjustable  arm  so  that  any  desired  rotor  sam- 
pling position,  YOj^,  could  be  obtained.  The  rotor  sampling  position 
was  measured,  in  degrees,  from  the  line  of  measurement  stations,  with 
the  stationary  blade  rows  set  as  indicated  earlier  (see  section  2.2.). 
Figure  5 shows  the  relative  position  of  the  stationary  blades  and 
several  rotor  sampling  positions.  The  sample-and-hold  circuit  was 
designed  and  constructed  by  the  Iowa  State  University  Engineering 

•k 

Research  Electronic  Shop  to  obtain  a 5 ysec  sample  with  each  revolu- 
tion of  the  rotor  shaft.  Two  scribe  marks,  one  on  the  rotating 
I portion  of  the  hub  surface  and  the  other  on  the  stationary  portion  of 

the  hub,  were  used  to  set  the  rotor  phase  lock  reference  position 
YOj^/Sj^  ” 0.0.  Once  the  reference  condition  was  established,  other 
rotor  sampling  positions  were  obtained,  by  moving  the  photoelectric 
pickup  with  respect  to  the  stationary  blade  .'iw  actuator  and  by 
observing  the  amount  of  movement  with  mechanical  scales.  The  reference 
setting  was  checked  daily  and  set  with  an  electronic  time  delay 
provision  in  the  sample-and-hold  circuit.  The  periodically  sampled 
signal  was  electronically  smoothed  with  the  signal  averaging  circuit, 
subsequently  read  by  the  digital  voltmeter,  and  finally  arithmetically 
averaged  and  stored  in  the  calculator  for  further  reduction  and  tape 
storage. 


* 

The  typical  rotor  wake  period  is  1.1  msec. 
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A desk-top  calculator  and  an  associated  multiple  channel 
scanning  digital  voltmeter  were  used  to  read,  store,  and  reduce  data. 

The  interfacing  between  the  calculator  and  voltmeter  enabled  the 
calculator  to  selectively  read  different  channels  and  store  readings. 

The  calculator  possessed  over  1400  memory  registers.  Additionally, 
with  tape  cassette  capability,  the  calculator  could  be  extended  to 
store  permanently  on  tape  any  data  and  programs  to  be  recalled  at  the 
user’s  convenience.  The  plotting  and  tabulating  of  results  were  done 
with  the  Iowa  State  University  Computation  Center  computing  system 
(IBM  360/65,  370/158). 


2.5.  Calibration  Nozzle 

An  air  nozzle  was  used  for  the  calibration  of  hot-wire  sensors. 

The  nozzle  has  a throat  diameter  of  25.4  mm  (1.0  in.)  and  a contraction 
ratio  of  144  to  1.  The  flow  at  the  nozzle  involved  a uniform  velocity 
profile  for  values  of  velocity  from  0.0  to  50  m/s.  Regulated  compressed 
air  provided  the  supply  air  and  a variable-current  heater,  blower,  and 
heat  exchanger  arrangement  was  used  for  air  temperature  control.  A 
special  probe  holder  permitted  the  sensing  portion  of  the  probe  to 
remain  in  the  same  position  in  the  nozzle  flow  while  the  probe  yaw  and 
pitch  angles  were  varied.  A telescope  was  used  to  visually  align  wires. 
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I 3.  EXPERIMENTAL  PROCEDURES  AND  DATA  REDUCTION 

I 

I All  measurements  presently  Involved  were  obtained  with  the  periodic- 

I average  measurement  system  described  earlier.  The  compressor  operating 

I 

I 

I point  was  held  constant  at  the  condition  indicated  in  Figure  7 

(1400±1  rpm,  0.42  flow  coefficient)  with  frequent  monitoring  and 
adjusting. 

The  standards  for  velocity,  pressure,  and  temperature  instru- 
ment calibrations  were  provided  by  the  calibration  nozzle,  the 
micromanoroeter , and  mercury- in- glass  thermometers.  Calibration  of 
all  electronic  equipment  was  performed  by  the  ISU  Engineering  Research 
Institute  Electronic  shop. 

The  desk-top  programmable  Hewlett-Packard  calculator  was  used 
extensively  for  calibration,  data  acquisition,  and  data  reduction. 
Programs  were  used  to  control  step-by-step  data  acquisition,  to  make 
preliminary  calculations,  to  print  out  data  and  preliminary  results, 
and  to  record  on  tapes  for  further  data  reduction.  Data  reduction 
programs  accepted  the  preliminary  results  and  produced  final  results. 

A list  and  brief  description  of  each  data  acquisition  and  reduction 
program  used  may  be  found  in  Appendix  A. 


3.1.  Periodic-Average  Sampling  Technique 

The  hot-wire  anemometer  output  signal  is  composed  of  a periodic 
component  and  a random  component.  With  controlled  sampling  and 
data  averaging  the  random  fluctuation  influence  can  be  made  as  small 
as  desired,  depending  on  the  nun^er  of  samples  taken  [4,7,18,23,25, 
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Figure  7.  Research  compressor  performance  curve  and 
operating  point. 

I 

r 


i' 


17 


26,28].  With  the  system  used  presently,  electronic  and  arithmetic 
averaging  of  data  could  be  accomplished.  The  electronic  averaging 
was  done  with  a low  pass  filter  having  a time  constant  of  1 sec. 
Schmidt  and  Okiishi  [28]  found  the  arithmetic  average  of  180  samples 
of  electronically-smoothed  data  to  be  approximately  equal  to  the 
arithmetic  average  of  1200  unfiltered  data  points.  More  current  data 
indicated  that  180  samples  of  elec Ironically- smoothed  data  were 
sufficient  to  make  the  random  component  of  the  unsteady  flow  negli- 
gibly small  except  at  the  IGV  row  exit  where  360  samples  were  needed 
because  of  the  separated  flow  activity  involved. 

The  data  sampling  system  was  triggered  by  a photoelectric  pickup, 
connected  to  the  circumferential  actuator  slider,  and  by  a single-slot 
disc  rotating  with  the  compressor  shaft.  The  photoelectric  pickup 
could  be  moved  relative  to  the  stationary  blade  rows  so  that  periodic- 
average  "snapshots"  of  the  flow  field  could  be  obtained  for  each  of 
several  time  sequenced  positions  of  the  sampled  rotor  blades.  Thus, 
periodic-average  flow  field  changes  could  be  measured  for  various 
sampled  rotor  positions.  This  kind  of  measurement  was  accomplished 
by  moving  the  stationary  blade  rows  and  the  photoelectric  pickup 
relative  to  the  stationary  probe  so  that  the  probe  was  effectively 
made  to  traverse  the  flow  field  over  one  stator  blade  pitch  spacing. 

3.2.  Single,  Slanted  Hot-Wire  Three-Dimensional 
Velocity  Measurement  Technique 

A single,  slanted  hot-wire  was  used  to  obtain  three-dimensional 
flow  field  parameters.  Before  the  measurement  technique  can  be 
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discussed,  some  relationships  linking  probe  geometry  and  hot-wire 
cooling  velocities  must  be  presented. 

3.2.1.  Probe  Geometry 

The  hot-wire  sensor,  the  probe  coordinate  system,  and  a general 
velocity  vector  are  shown  in  Figure  8.  The  coordinate  system  is 
fixed  to  the  probe  with  the  x-z  plane  lying  on  the  sensing  portion 
of  the  probe  and  the  probe  axis  and  with  the  y-axis  perpendicular  to 


the  x-z  plane  and  centered  on  the  sensor.  The  wire  is  slanted  at  an 

angle,  9^,  to  the  x-axis.  The  velocity  vector,  V,  can  be  resolved 

into  orthogonal  components  along  x,  y,  and  z for  each  orientation  of 

the  wire.  When  the  probe  and  its  coordinate  system  are  rotated 

about  the  z-axis,  the  projected  yaw  angle,  6^,  changes  by  the  amount 

of  turning,  whereas,  the  pitch,  9 , remains  the  same.  The  sensor 

P 

yaw  angle,  a,  was  defined  as  the  angle  between  the  intersection  of 

the  probe  sensor  unit  vector,  X,  and  the  velocity  vector,  V.  To 

obtain  a relationship  between  a and  9_,  9 , and  9 , the  dot  product 

0 p y 

of  the  two  vectors  is  taken: 


A = cos  9q  i + sin  9q  k 


V = -V  cos  9 cos  9 X - V cos  9 sin  9 1 - V sin  9 k (2) 

p y p y p 

A*V  = IaI  lv|  cos(180  - nt)  = -'vl  cos  9.  cos  9 cos  9 
' 0 p y 

-Ivl  sin  9 sin  9 (3) 

U p 

/.cos  a * cos  cos  0 cos  0 + sin  0^  sin  0 (4) 

0 p y Op 
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Figure  8.  Hot-wire  configuration  relating  velocity  vector,  V,  to 
hot-wire  sensor  and  probe  coordinates  x,  y,  z. 
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3.2.2.  Effective  Cooling  Velocity  Ratio 

The  hot-wire  anemometer  output  was  correlated  to  fluid  velocity 
with  a fourth-order  fit  so  that  the  anemometer  signal  could  elec- 
tronically be  linearized.  This  linearized  signal,  E^,  was  correlated 
to  velocity  with  the  sensor  angle  normal  to  flow  to  obtain 

V . Kj  + KjE,  + (5.) 

where  V is  the  absolute  fluid  velocity.  The  constants  K^,  K^,  and 
were  determined  with  a least  squares  fit.  Whenever  the  probe  is 
oriented  to  the  flow  at  other  than  a sensor  yaw  angle,  a,  of  90°,  the 
velocity  read  can  be  considered  an  effective  cooling  velocity,  V^, 
where 

. Ki  * ,5^) 

The  presently  used  velocity  measurement  technique  was  based  on  knowing 
a precise  relationship  for  the  effective  cooling  velocity/absolute 
velocity  ratio,  V^/V,  for  various  orientations  of  the  probe  in  the 
flow  stream. 

Experiments  conducted  by  Schmidt  and  Okiishi  [28]  showed  this 
velocity  ratio  was  strongly  dependent  on  sensor  yaw  angle  and  weakly 
dependent  on  pitch  angle  and  velocity  level.  The  recommended  correlation 
is  as  follows: 

V /V  = b_  + b.a  + b„9  + b_V  + b.a^  + b,6  + b,V 

e 01  2p3  4 5p  o 

+ b.,ae  + b-aV  + b-0  V (6) 

/ p o 9 p 

The  coefficients  b^  through  b^  were  determined  from  a least  squares 
fit  of  data  as  described  in  the  calibration  procedure  section  (see 
section  3.3.3.). 


3' 

V 
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3.2.3.  Measurement  Technique 

The  wire  was  rotated  in  the  flow  field  to  three  different  positions 
denoted  as  a,  b,  and  c (see  Figure  9).  These  probe  positions  relate  to 


yaw  angles  of  9 ,6  , , and  0 which  were  set  as  indicated  below: 

y,a  y,b  y,c 


9 


9 


y.a  y 

6 , = 0 - m, 

y.b  y b 


0 = 9 - m 

y.c  y c 


(7) 

(8) 
(9) 


where  and  are  probe  turning  angle  increments  from  the  a location. 
For  each  orientation  of  the  wire,  effectively  1200  samples  were  taken 
and  electronically  and  arithmetically  averaged  to  obtain  a time-average 
effective  cooling  velocity  for  that  particular  probe  position.  For 
each  periodic-average  velocity  vector  measurement,  a total  of  six 
equations  like  Eqs.  (4)  and  (6)  were  obtained  for  the  orientations  of 
the  wire.  These  equations  are: 

For  position  a 

V /V  = b„  + b,a  + b^0  + b,V  + b,a  ^ + b^0  ^ 

e,a  0 la  2p  3 4a  5p 


+ b,V  + b,a  0 + b-a  V + b-.9  V 

6 7ap  8a  9p 

cos  a « cos  9-  cos  6 cos  6 + sin  9„  sin  6 

a Op  y,a  0 p 

For  position  b 

■ ”0  “1%  * ''2«p  * '>3''  ^ 

+ b,\^  + b.,a,0  + b„a.V  + b„0  V 


(10) 

(11) 


’7“b"p  ■’8%' 


9 p' 


COS 


“b 


cos  0_  cos  0 cos  0 . + sin  0.  sin  0 

0 p y,b  Op 


(12) 

(13) 
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For  position  c 

V /V  = b-  + b.a  + b_0  + b,V  + b,a  ^ + b^0  ^ 

e,c  0 Ic  2p  3 4c  5p 

+ b,V^  + b_a  e + b_a  V + b-0  V 
6 7cp  8c  9p 

cos  a = cos  0„  cos  0 cos  0 + sin  0^  sin  0 

c Op  y,c  0 p 


(14) 

(15) 


By  substituting  Eqs.  (7),  (8),  and  (9)  into  Eqs.  (11),  (13),  and  (15), 

the  six  unknown  variables  a , a.  , a , 0 , 0 , and  V remain  in  the  six 

a b c p y 

equations,  (10)  through  (15).  These  equations  were  solved  simultaneously 
using  the  Newton-Raphson  method.  The  three-dimensional  velocity  vector 
is  completely  described  with  the  variables  0,0,  and  V known 

p y 

relative  to  the  probe  position,  a. 

To  set  the  three  wire  orientations,  a minimum  effective  velocity 

angle,  6 , was  obtained  as  shown  in  Fiure  9.  From  this  minimum 

mv 

effective  velocity  angle  position,  the  probe  was  rotated  to  positions 
a,  b,  and  c according  to  the  following  schedule: 

0 = 20° 

a,  off 

0K  = 60° 

b, off 

0 _ = -20° 

c, off 

From  Figure  9,  it  can  be  seen  that 


0 

0 

= -20° 

y.a 

y 

0 . * 
y.b 

<I> 

- m^  ==  -60' 

0 

0 

- m =20° 

y.c 

y 

c 

40'- 


m = -40 
c 


and 
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This  positioning  of  the  probe  was  done  for  each  velocity  vector 
obtained.  For  all  of  the  measurements,  the  probe  was  positioned  so 
that  the  longer  wire  support  was  downstream  of  the  sensing  portion 
of  the  wire  to  minimize  prong  interference. 


3.3.  Calibration  Procedures 


The  calibration  nozzle  was  used  to  provide  a uniform  velocity  air 
stream  for  velocity  calibration.  Plenum  pressures  were  observed  with 
Inclined  water-in-glass  manometers.  Velocity  was  calculated  using  the 
following  equation: 


where 


V 


AP 

n 


(16) 


V » velocity,  m/s 

2 

g^  = gravitational  constant  1.0  kg  m/Ns 

3 

Yjj  Q = specific  weight  of  water,  N/m 

3 

0 = density  of  air,  kg/ro 

AP  = differential  pressure  between  plenum  pressure  and 
atmospheric  pressure,  meters  of  water 

Calibration  programs  were  used  with  the  desk-top  calculator. 

Sufficient  warm-up  time  was  allowed  to  obtain  steady  state  conditions. 

Fluid  temperatures,  barometric  conditions,  and  test  conditions  were 

entered  into  the  calculator  before  and  after  each  calibration.  Three 

hot-wire  calibrations  were  made:  linearlzer  velocity  calibration, 

second  order  velocity  calibration,  and  effective  cooling  velocity  ratio 


calibration. 
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3.3.1.  Linearlzer  Velocity  Calibration 

The  linearlzer  was  used  to  approximate  the  anemometer  output  with 
a fourth  order  polynomial.  The  "zero"  degree  term  of  the  polynomial 
was  set  equal  to  zero,  and  coefficients  were  determined  by  a least 
squares  curve  fit,  thus  producing  a linear  relationship  between  velocity 
and  anemometer  output  with  the  wire  positioned  normal  to  the  flow  of 
the  calibration  nozzle. 

The  wire  was  optically  positioned  6.35  mm  (0.25  in.)  away  from 

the  nozzle  exit  plane.  Wit.i  the  air  velocity  set  close  to  the  average 

value  of  the  range  of  velocities  anticipated,  the  cold  combined 

resistance  of  the  sensor,  cables,  and  leads  was  measured,  R The 

s,c,d 

operating  resistance  of  the  sensor  was  calculated  with  an  overheat 

ratio  of  1.8.  Because  the  cablt  and  lead  resistances  were  not  involved 

with  the  overheat  calculation,  they  were  subtracted  from  the  cold 

resistance  reading  at  the  anemometer  deck  and  then  subsequently  added 

back  on  as  demonstrated  in  Eq.  (17). 

R ,=1.8(R  .-R.-P.-RJ 

s,op,d  s,c,d  cb  ph  pi 

+R,+R^+R,  (17) 

cb  ph  pi 

where 


' op  d ~ setting  of  sensor  operating  resistance,  ohn 

_ ^ ^ = deck  reading  of  combined  sensor,  cable,  and  lead 
’’  ’ resistance,  ohms 


R , = cable  resistance,  ohms 

CD 

Rpj^  = probe  holder  resistance,  ohms 

R , = probe  lead  resistance,  ohms 
pi 


r 
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Probe  holder  and  cable  resistance  were  measured  with  an  impedance 
bridge  while  the  probe  lead  resistance  was  taken  as  furnished  by  the 
manufacturer.  The  linearizer  calibration  system  included  the  hot-wire, 
anemometer,  DVM,  and  desk-top  calculator  only.  The  linearizer  was 
correlated  over  a range  of  nozzle  air  velocities  (0.0  to  23  m/s)  in 
fifteen  increments.  The  linearized  polynomial  velocity  was  compared 
to  the  actual  nozzle  velocity.  The  acceptance  criterion  for  the  points 
was  that  the  error  should  be  less  than  or  equal  to  1%.  Linearizer 
calibration  was  performed  with  every  effective  cooling  velocity  ratio 
calibration. 

3.3.2.  Second  Order  Velocity  Calibration 

A second  order  correlation,  Eq.  (5),  between  the  linearized 
anemometer  output  and  velocity  level  with  the  wire  normal  to  the  flow 
was  performed  before  and  during  effective  cooling  velocity  calibration 
and  data  acquisition  to  avoid  error  due  to  electronic  equipment  drift. 

The  second  order  calibration  was  done  with  the  wire  positioned 
normal  to  the  nozzle  flow.  Any  other  wire  yaw  angle  does  not  indicate 
an  absolute  velocity  but  rather  an  effective  cooling  velocity,  V^. 

The  second  order  calibration  system  included  the  probe,  anemometer, 
linearizer,  DVM,  and  desk-top  calculator.  Fourteen  points  were 
taken  for  each  calibration  for  air  velocities  ranging  from  4 to  23  m/s. 
Acceptable  error  was  confined  to  less  than  1%  in  the  velocity  working 
range . 

3.3.3.  Effective  Cooling  Velocity  Calibration 

The  effective  cooling  velocity  calibration  was  performed  to 
obtain  the  coefficients  in  the  effective  cooling  velocity  ratio  equation: 


and  pitch  angles  and  velocities. 

The  probe  was  positioned  in  the  nozzle  jet  stream  in  the  same 
manner  as  for  the  linearizer  calibration.  For  a particular  velocity 
level,  the  probe  was  set  at  each  of  six  pitch  angles  and  rotated 
from  0 to  90°  in  yaw  angle  increments  of  5°.  This  procedure  was 
continued  until  data  for  all  of  the  velocity  levels  (11.6,  15.2, 

19.2,  22.3  m/s)  and  pitch  angles  (-9  to  6°,  in  increments  of  3°)  were 
acquired.  This  procedure  was  also  used  for  probe  yaw  angles  of  0 to 
-90°  because  of  the  lack  of  perfect  symmetry  wire  response  to  yaw 
angle.  Thus,  during  data  reduction,  two  sets  of  coefficients  (b^  through 
bg)  were  used  depending  on  the  yaw  orientation  of  the  wire.  Positions 
a and  b involved  the  0 to  -90  calibration  coefficients  while  position 
c used  the  0 to  +90  calibration  coefficients. 

After  all  the  data  were  taken,  the  coefficients  were  solved  using 
a least  squares  method.  Acceptable  error  was  limited  to  being  less 
than  2%  with  the  majority  of  the  errors  less  than  1%. 

3.4.  Data  Acquisition 

The  data  acquisition  system  used  was  presented  earlier  in  schematic 
form  in  Figure  6.  Prior  to  taking  data,  equipment  and  room  temperature 
equilibrium  were  established,  inclined  manometers  zeroed,  cold  resis- 


J 


tance  of  the  wire  at  the  deck  set,  flow  coefficient  set,  linearizer 


28 


coefficients  adjusted,  second  order  calibration  performed,  probe 

actuator  positioned  in  the  compressor,  and  radial  location  set.  After 

setting  and  adjusting  the  time  delay  device  for  YO_/S  =0.0,  the 

R R 

desired  rotor  sampling  position  could  be  set  by  adjusting  the  photo- 
electric pickup  position  relative  to  the  stationary  blade  rows. 

With  the  probe  sensor  positioned  at  90°  (normal  to  the  compressor 
axis),  a trace  of  the  blade-to-blade  variation  of  effective  cooling 
velocity  was  made  with  an  x-y  storage  oscilloscope  where  x was  the 
potentiometer  reading  from  the  circumferential  actuator,  and  y was 
the  time-averaged  hot-wire  sampled  signal.  This  survey  was  made  by 
moving  the  stationary  blade  rows  and  photoelectric  pickup  relative 
to  the  probe.  By  using  this  procedure,  the  probe  was  effectively  mide 
to  traverse  the  flow  field  for  one  stator  blade  pitch.  From  this 
x-y  trace,  the  position  of  the  wake  could  be  detected  and  test  points 
appropriately  scheduled.  The  circumferential  actuator  system  was 
then  set  back  to  zero,  Y/S^  = 0.0,  for  the  first  test  data  point.  With 
another  x-y  storage  oscilloscope,  a trace  of  effective  cooling  velocity 
against  probe  yaw  angle  was  repeatedly  obtained  by  rotating  the  probe 
about  its  axis.  In  this  fashion,  a minimum  effective  cooling  velocity 
angle,  6^^,  could  be  discerned.  From  this  minimum  angle  position,  the 
probe  could  be  appropriately  positioned  to  locations  a,  b,  and  c, 
discussed  earlier  (see  Figure  9).  After  periodically  sampled  data 
wore  obtained,  with  the  probe  at  each  position  (a,  b,  and  c),  the 
circumferential  actuator  was  shifted  to  a new  position,  Y/S^,  and  the 
above  described  procedure  was  repeated.  Thirty  velocity  vectors  from 
blade-to-blade  were  measured  for  each  test  run.  After  the  last  velocity 
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vector  of  a set  of  thirty  was  measured,  the  acquired  data  were  stored 
on  tape  for  further  reduction.  The  probe  was  then  removed  from  the 
compressor  and  a new  second  order  velocity  calibration  was  made.  Second 
order  velocity  calibrations  were  performed  before  each  set  of  thirty 
velocity  vectors  was  measured. 


3.5.  Data  Reduction 


Data  reduction  programs  for  the  desk-top  calculator  were  used  to 

accept  and  reduce  stored  data.  The  six  equations,  Eqs.  (10)  through 

(15),  were  solved  simultaneously  with  the  Newton-Raphson  method. 

Usually,  less  than  five  iterations  were  necessary  for  convergence. 

The  compressor  coordinate  system  is  shown  in  Figure  10.  The  R 

coordinate  direction  is  positive  outward  from  the  hub,  while  the  9 

coordinate  is  positive  in  the  direction  of  rotation.  The  Z-axis  is 

aligned  axially,  positive  in  the  direction  of  fluid  motion.  Also 

shown  are  the  sign  conventions  relating  to  6,,  6 , V , V , and  V„. 

8 r r z 0 

With  the  probe  immersed  in  the  compressor  annulus,  the  probe  coordinates 
X and  y (see  Figure  8)  are  always  in  the  same  plane  as  Y and  Z of 
the  compressor  coordinate  system,  and  z from  the  probe  coordinate 
system  coincides  with  the  R compressor  coordinate  direction. 

The  calculated  pitch  angle  and  radial  angle  are  related  by: 


6,  - -e 

r P 


The  tangential  angle  is  related  to  the  calculated  yaw,  0^,  by: 

6o  - 6 +9  + e 

0 mv  a, off  y 


three-dimensional  periodic-average  velocity  and  angle  parameters 
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The  values  of  S and  9 ..  were  known  from  data  acquisition, 

mv  a , o f f ^ 

other  flow  parameters  calculated  are  presented  below: 

(1)  Axial  velocity,  m/s,  Eq.  (B-15) 

(2)  Absolute  tangential  velocity,  m/s,  Eq.  (B-16) 

(3)  Radial  velocity,  m/s,  Eq.  (B-14) 

(A)  Relative  velocity,  m/s,  Eq.  (B-17) 

(5)  Relative  tangential  velocity,  m/s,  Eq.  (B-18) 

(6)  Relative  tangential  angle,  degrees,  Eq . (B-19) 

A complete  list  of  equations  used  appears  in  Appendix  B. 
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4.  PRESENTATION  AND  DISCUSSION  OF  DATA 


The  test  results  are  presented  and  discussed  in  this  section. 
Primary  flow  variables  involved  are  tabulated  in  Appendix  C.  Velocity 
and  flow  angle  component  variation  graphs  (scalar  plots)  are  used  to 
present  the  periodic-average  data  (see  Figures  11,  12,  and  13).  The 
scalar  plots  include  axial,  tangential,  and  radial  velocities  supple- 
mented with  tangential,  radial,  and,  when  useful,  relative  tangential 
angles.  Uncertainty  and  accuracy  are  previously  discussed  by  Schmidt 
and  Oklishi  [28].  To  reflect  and  facilitate  data  interpretation, 
blade- to-blade  plane  flow  plots  (cascade  plots),  s.iowing  the  periodic- 
average  location  of  blade  wake  flows  at  different  radii  for  several 
distinct  orientations  of  the  rotor  blades,  were  constructed  (see 
Figure  14).  The  procedure  used  to  construct  the  cascade  plots  is 
explained  in  detail  below.  A physical  explanation  of  the  wake  trans- 
port and  interaction  effects  involved  is  proposed. 

4.1.  Construction  of  Cascade  Wake  Plots 


Statiuiiary  blade  section  profiles  were  positioned  on  the  plots  to 
scale  according  to  the  minimum  sound  schedule  (see  section  2.2.). 

The  rotor  blade  sections  were  located  at  their  periodic  sampling 
positions. 

The  obvious  wake  region  locations  suggested  by  the  periodic-average 
flow  data  were  drawn  on  the  cascade  plots  first.  It  should  be  noted 
that  periodic-average  flow  data  collection  was  limited  to  thirty  points 
along  the  broken  lines  connecting  crosses  between  blade  rows,  Y/S  = 
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(a)  (Continued). 
Figure  12.  (Continued) 
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Figure  12.  (Continued) 
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Figure  12.  (Continued) 
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lOZ  passage  height  from  hub.  30%  passage  height  from  hub. 
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(h)  First  stator  exit  flow  at  rotor  sampling  position  YO_/S„  - 0.69  and  five  different 
passage  heights. 
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0.0  to  1.0  (see  Figure  14).  Behind  the  IGV  blade  row,  the  IGV  wake 
regions  were  easily  identified  and  located  on  the  scalar  plots  for 
all  radial  span  locations  and  rotor  sampling  positions.  Behind  the 
first  rotor,  the  rotor  wake  regions  were  easily  identified  at  all  radii 
and  rotor  sampling  positions,  but  the  IGV  wake  regions  were  more 
difficult  to  see  and  were  discernible  at  some  span  locations  only. 

Behind  the  first  stator,  the  stator  wake  regions  were  clearly  identi- 
fiable, but  the  rotor  wake  and  IGV  wake  effects  were  not  clearly 
distinguishable  in  all  cases.  At  the  stator  exit,  the  rotor  wake 
region  could  often  be  identified  by  noting  the  characteristic  predominant- 
ly-outward  radial  flow  Involved,  for  example,  as  seen  in  Figure  13(a) 
for  YOj^/Sj^  = 0.0  near  Y/Sg  = 0.2.  Rotor  wake  effects  behind  the  stator 
row  were  connected  to  corresponding  rotor  blades  consistently  with 
the  time-average  data  reported  in  Reference  28.  In  order  to  ascertain 
which  IGV  blades  caused  the  IGV  wake  effects  observed  behind  the  first 
rotor  and  the  first  stator,  an  experiment,  inspired  by  Smith's  [29]  dye 
source  analogy,  was  conducted  as  described  below.  The  sensing  of 
heated  air,  flowing  in  the  wake  of  a pen-type  soldering  iron  held 
upstream  of  the  IGV  row,  with  a thermocouple  immersed  downstream  of 
the  first  rotor  row  and  the  first  stator  row  revealed  the  match  between 
downstream  wake  effect  with  upstream  wake  production.  It  was  assumed 
that  in  a time-average  sense,  the  heated  soldering  iron  wake  flow 
behaved  like  the  IGV  wake  flow,  and,  thus,  the  IGV  wake  "avenues", 
wliich  can  be  described  as  stationary  paths  within  which  the  chopped 
wakes  move  downstream,  were  conceived  to  be  parallel  to  the  heated  air 
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flow  paths  (see  Figure  15).  The  procedure  described  above  resulted 
in  cascade  plots  in  skeletal  form  only.  Completion  of  the  cascade  plots 
involved  careful  interpolation  between  cascade  wake  plots  for  the  same 
radii  and  different  rotor  sampling  positions  and  cascade  plots  for  the 
same  rotor  sampling  position  but  different  radii.  By  using  a wake 
template  involving  an  increasing  wake  width,  the  wake  thickening  effect, 
suggested  in  References  2 and  3,  was  approximated.  The  chopped 
segments  were  caused  to  reflect  Smith's  [29]  chopping  effects  by  rotat- 
ing each  bounded  wake  segment  (bounded  either  by  blade  surfaces  or 
wake  regions)  around  its  center.  The  usefulness  of  these  cascade  plots 
should  become  evident  as  the  test  results  are  discussed  in  the  next 
section. 

4.2.  Discussion  of  Data 

The  results  of  the  periodic-average  measurements  made  behind  the 
IGV  row  are  presented  in  velocity  and  angle  component  form  in  Figure 
11  and  will  be  discussed  first.  The  two  plots  of  axial  velocity 
variation  from  blade-to-blade  at  50%  passage  height  from  the  hub  (PHH) , 
for  two  different  rotor  sampling  positions,  indicate  that  the  rotor 
potential  flow  field  could  influence  the  flow  at  the  IGV  row  exit 
measurement  plane.  Walker  and  Oliver  [3]  suggest  that  3%  to  4% 
potential  flow  effect  can  be  expected  for  this  particular  axial  spacing 
between  blade  rows.  The  present  data  are  in  close  agreement  with  this 
observation.  Rotor  leading  edge  locations,  as  seen  in  Figure  16  for 
rotor  sampling  positions,  YO^^/Sj^  = 0.0  and  0.34,  relate  to  the  axial 
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Figure  15.  Time-average  heated  air  avenue  location  in  cascade. 
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(a)  Rotor  sampling  position  YO^^/SR  •0.0. 
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(b)  Rotor  sampling  position  YO^/S^  ■ 0.34. 

Figure  16.  Relative  locations  of  IGV  and  rotor  blades  for  two 
rotor  sampling  positions  at  502  span. 
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velocity  deficit  locations,  Y/S„  = 0.5  for  Y0„/S„  = 0.0,  and  Y/S„  =0.3 

S K K o 

for  YOj^/Sj^  = 0.34.  The  rotor  upstream  effect  is  more  apparent  in 

the  tangential  velocity  and  angle  variations  from  blade-to-blade  for 

rotor  san^jling  position,  YO  /S_  = 0.34.  Radial  velocities  and  angles 

do  not  reveal  any  rotor  upstream  effect  when  the  rotor  leading  edge 

is  located  in  the  freestream  portion  of  the  IGV  exit  flow  (Y0^/S„  = 

0.34).  However,  a general  increase  in  outward  radial  flow  of  the  IGV 

wake  fluid  is  evident  when  the  rotor  leading  edge  is  downstream  from 

the  IGV  trailing  edge  (Y0_/S  = 0.0).  Examination  of  the  IGV  exit 

K R 

flow  fields  from  hub  to  tip  for  rotor  sampling  position,  YO  /s  = 0.0, 
only,  suggests  that  these  wakes  becoti^  narrower  from  10%  to  90%  PHH. 

At  70%  span,  the  IGV  wake  is  curiously  only  half  as  deep  as  the  wakes 
at  50%  and  90%  PHH.  There  is  a general  decrease  in  freestream  axial 
velocity  for  90%  PHH,  indicating  an  outer  wall  boundary  layer  effect. 

It  has  been  observed  that  for  rotor  sampling  position,  YO  /S_  = 0.0, 
the  rotor  leading  edge  is  directly  downstream  from  the  IGV  trailing 
edge.  Also,  the  rotor  leading  edge  is  closer,  axially,  to  the  IGV  exit 
measurement  plane  at  10%  PHH  and  is  farther  at  90%  PHH  (see  Figure  4) . 
Tangential  velocities  follow  a predictable  and  systematic  pattern  of 
smaller  tangential  velocities  in  the  IGV  wake.  Tangential  angles,  on 
the  other  hand,  generally  tend  to  decrease,  suddenly  Increase,  and  then 
decrease  again  in  proceeding  from  the  suction  side  to  the  pressure  side 
of  the  IGV  wake.  Radial  velocities  and  angles  generally  exhibit  a 
curious  outward  flow  in  the  wake,  except  at  50%  and  70%  PHH.  Irregular 
separated  flow  activity  was  observed  on  the  suction  side  of  the  IGV 
wake  during  data  acquisition.  There  did  not  appear  to  be  any  difference 
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In  separated  flow  Irregularity  from  hub  to  tip.  To  obtain  meaningful 
periodic-average  data,  approximately  2400  Instead  of  the  usual  1200, 
effective  samples  were  obtained  for  each  orientation  of  the  wire. 

First  rotor  periodic-average  exit  flow  parameters  for  50%  span 
are  presented  in  Figure  12(a)  for  six  rotor  sampling  positions  and 
will  be  discussed  next.  Corresponding  cascade  plots  are  displayed 
in  Figure  14(c).  Note  that  the  rotor  wake  regions  shown  in  the  sequence 
of  Figure  12 (a)  are  not  from  the  same  rotor  blade  but  rather  are 
from  two  adjacent  rotor  blades.  This  fact  may  be  verified  by  following 
the  time-sequential  movement  of  the  rotor  blades  in  the  cascade  plots 
(Figure  14 [c]).  The  chopped  IGV  wakes  seen  in  the  cascade  plots 
follow  a fixed  avenue.  Within  this  avenue,  the  wake  segments  move 
downstream,  as  indicated  in  Figure  14(c).  Periodically,  the  rotor  wake 
region  can  occupy  the  same  portion  of  the  measuring  plane  between 
Y/S  = 0.0  and  1.0  as  the  IGV  wake  avenue.  The  deepest  rotor  wake 
region  measured  was  for  a rotor  sampling  position  of  Y0^/S„  = 0.69. 

For  this  rotor  position,  the  cascade  plot  (Figure  14 [c])  indicates 
that  the  rotor  wake/lGV  wake  interaction  occurs  slightly  upstream 
from  the  measurement  plane.  For  this  rotor  position,  the  circumferential 
variation  of  flow  parameters  in  the  rotor  wake  region  is  characterized 
by  larger  deviation  and  absolute  flow  angles  and  by  smaller  absolute 
tangential  velocities  and  axial  velocities  than  observed  in  the  rotor 
wake  regions  associated  with  other  rotor  sampling  positions  where  the 
rotor  wake/IGV  wake  interaction  occurs  closer  to,  at,  or  downstream 
of  the  measurement  plane.  When  the  rotor  wake/lGV  wake  Interaction 
occurs  slightly  upstream  or  at  the  measurement  plane,  the  wake  region 


will  hereafter  be  called  an  "interacted  wake  region."  When  the  rotor 

wake/lGV  wake  interaction  occurs  downstream  of  the  measurement  plane, 

the  wake  region  will  hereafter  be  referred  to  as  a "noninteracting 

wake  region."  The  dependence  of  rotor  wake  region  behavior  on  rotor 

sampling  position  is  explainable  in  terms  of  the  physical  reasoning 

proposed  by  Kerrebrock  and  Mikolajczak  [2]  which  involved  the 

relative  flow  of  chopped  stator  wake  fluid  moving  toward  the  rotor 

pressure  surface  due  to  slip  (see  Figure  17[a]).  Such  flow  in  a 

chopped  IGV  wake  would  tend  to  result  in  locally  reduced  axial  velocities 

and  Increased  deviation  angles  on  the  rotor  suction  surface  side  and 

in  an  accumulation  of  more  slower-moving  fluid  on  the  rotor  pressure 

surface  side.  The  mechanism  for  these  interaction  effects  is  best 

explained  in  this  way:  The  collected  IGV  wake  fluid  "builds  up" 

a surplus  of  fluid  on  the  pressure  side  of  the  rotor  blade,  while  on 

the  suction  side  of  the  rotor  blade,  fluid  is  drawn  towards  the  rotor 

pressure  surface  of  the  adjacent  blade.  Thus,  a narrower  wake  region, 

larger  deviation  angles,  and  smaller  axial  velocities  occur  at  the 

measurement  plane  when  the  suction  side  wake  interaction  effects 

predominate  there,  as  for  YO  /S  • 0.69.  For  YO  /S  » 0.69,  the 

R R R R 

smaller  absolute  tangential  velocities  and  larger  absolute  tangential 

angles  are  explained  with  the  help  of  a velocity  vector  diagram 

(Figure  18),  where  three  velocity  triangles  are  drawn,  one  each  for 

free  stream,  noninteracting  wake,  and  interacting  wake  regions.  For 

rotor  sampling  position  YO  /S„  - 0.83,  the  rotor  blade  pressure 

surface  side  interaction  effect  is  stronger  at  the  measurement  plane, 

resulting  in  a wake  zone  spreading  effect.  For  YO  /S  ■ 0.0,  both 

R R 


Figure  18.  Plane  velocity  vector  triangles  for  fluid  in  an  interacted 
wake,  noninteracted  wake,  and  free  stream  for  the  first 
rotor  exit  flow  at  midspan. 
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pressure  and  suction  surface  interaction  effects  are  noticeable  at 
the  measurement  plane,  the  result  being  a broad  shallow  wake  region. 
Radial  velocity  and  angle  distribution  change  trends  for  different 
rotor  sampling  positions  are  not  apparent  from  the  data. 

Prior  to  measuring  the  details  of  the  flow  behind  the  first  rotor 
at  other  radial  locations  in  the  compressor  annulus,  approximate 
periodic-average  hot-wire  oscilloscope  traces  of  a combination  of 
axial  and  radial  velocities  were  studied  to  determine  which  rotor 
sampling  positions  should  be  selected  to  show  the  largest  changes  in 
periodic-average  flow  due  to  rotor  wake/lGV  wake  interaction.  For 
radial  locations  10%  and  90%  PHH,  the  appropriate  rotor  sampling 
positions  to  use  were  not  obvious  because  the  approximate  data  were 
not  definitive  enough.  At  the  other  radial  positions,  more  reasonable 
rotor  sampling  position  selections  could  be  made  with  the  approximate 
data. 

The  rotor  exit  data  for  10%  PHH  for  four  rotor  sampling  positions 
are  shown  in  Figure  12(b).  The  related  cascade  plots  are  displayed 
in  Figure  14(a).  Unfortunately,  the  data  involved  only  rotor  wake/ 

IGV  wake  interaction  at  the  measurement  plane  for  rotor  sampling 
positions,  YO^/S^  = 0.0,  0.32,  0.34,  and  downstream  of  the  measurement 
plane  for  YO^/Sj^  = 0.69.  Data  for  interaction  occurring  slightly 
upstream  of  the  measurement  plane,  which  might  have  Indicated  a com- 
paratively deeper  rotor  wake  region,  were  not  obtained.  The  information 
shown  is  consistent  with  data  already  discussed. 

The  periodic-average  data  for  30%  PHH  are  shown  in  Figure  12(c) 


for  rotor  sampling  positions  YOj^/Sj^  =0.0  and  0.42,  which  represent 
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interacted  and  noninteracted  wake  regions  respectively.  The  cascade 
plots  (Figure  14[b])  indicate  that,  for  rotor  sampling  position 
YO  /S  = 0.0,  a characteristically  deeper  interacted  wake  region  is 
produced  with  the  rotor  wake/lGV  wake  interaction  occurring  upstream 
of  the  measurement  plane.  The  flow  parameters  exhibit  the  same  charac- 
teristics as  did  the  rotor  exit  flow  for  50%  span,  Y0_/S„  = 0.69, 

R R 

namely,  reduced  axial  velocities  and  absolute  tangential  velocities 
and  increased  deviation  angles  and  absolute  tangential  angles  in  the 
interacted  wake  region. 

First  rotor  exit  flow  component  plots  for  70%  PHH  and  two  rotor 

sampling  positions,  YO^^/S^^  = 0.0  and  0.53,  are  shown  ir  Figure  12(d). 

The  corresponding  cascade  plots  (Figure  14[d])  reflect  an  interacted 

rotor  wake  (interaction  occurring  upstream  of  the  measurement  plane) 

for  YOj^/Sj^  = 0.53  and  a noninteracted  wake  for  rotor  sampling  position 

YOj^/Sj^  = 0.0.  For  YO^/S^  = 0.53,  the  rotor  wake/IGV  wake  interaction 

effects  are  typical  of  a deeper  interacted  rotor  wake  region  except  for 

the  tangential  angles  which  are  slightly  smaller  in  the  interacted 

wake  region  than  in  the  noninteracted  region.  This  inconsistency  is 

probably  due  to  the  fact  that  the  interacted  wake  region  for  YO  /S  = 

R R 

0.53  is  not  the  deepest  possible  and,  thus,  does  not  totally  possess  all 
of  the  characteristics  exhibited  by  such  wake  regions. 

The  data  for  90%  PHH  rotor  exit  flow  at  two  rotor  sampling 
positions,  YO^^/Sj^  “0.0  and  0.63,  are  shown  in  Figures  12(e)  and  14(e). 
The  rotor  wake/IGV  wake  Interaction  occurs  at  the  measurement  plane  for 
YOj^/Sj^  - 0.63  and  slightly  downstream  for  YO^/S^  » 0.0,  The  rotor 
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wake  region  corresponding  to  interaction  upstream  of  the  measurement 
plane  was  not  observed.  The  trends  indicated  by  the  present  data 
are  consistent  with  similar  measurements  made  elsewhere  in  the 
compressor  annulus. 

Rotor  exit  flow  data  for  seven  radial  positions  but  only  one  rotor 

sampling  position,  YO  /S  = 0.0,  are  presented  in  Figure  12(f). 

K K 

Axial  velocity  variations  indicate  a general  decrease  in  wake  depth 
from  hub  to  tip  and  a circumferential  shifting  consistent  with  change 
of  blade  twist  (see  Figure  19)  except  at  95%.  The  curious  circum- 
ferential shift  at  95%  PHH  is  not  clearly  explainable,  but  it  possibly 
is  due  to  tip  scraping  and  clearance  leakage  effects.  Radial  velo- 
cities and  angles  show  tip  scraping  and  clearance  leakage  effects  at 
90%  and  95%  PHH.  The  outer  casing  boundary  layer  extends  to  95%  PHH. 

A hub  boundary  layer  is  not  apparent  in  the  velocity  data  for  5%  PHH 
although  radial  velocities  and  angles  do  indicate  the  presence  of 
secondary  flows.  At  10%  and  30%  PHH,  deep  interacted  wake  regions  are 
observed  while  shallower  noninteracted  wake  regions  are  seen  at  70% 
and  90%  PHH.  IGV  wake  effects  seem  to  be  more  apparent  toward  the 
tip,  possibly  due  to  the  smaller  amount  of  rotation  of  the  chopped  IGV 
wakes  involved  there. 

First  stator  periodic-average  exit  flow  data  for  50%  PHH  and  six 
rotor  sampling  positions  are  presented  in  Figure  13(a)  with  associated 
cascade  plots  in  Figure  14(c).  It  is  obvious  that  the  first  stator 


exit  flow  is  much  more  complicated  than  the  first  rotor  exit  flow 
because  of  the  stator,  rotor,  and  IGV  wake  interactions  involved.  The 
stator  exit  flow  variation  with  rotor  sampling  position  will  be  discussed 


90%  PHH 

50%  PHH 

10%  PHH 


-h 

Figure  19.  Compressor  rotor  blade  sections  at  hub,  midspan 
and  tip  locations. 
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first.  As  mentioned  by  Kerrebrock  and  Mikolajczak  [2],  the  chopped 

rotor  wake  fluid  tends  to  move  from  the  stator  suction  surface  to 

the  pressure  surface  because  of  the  slip  motion  in  the  wake  (see 

Figure  17 [b]).  Such  flow  in  a chopped  rotor  wake  would  tend  to 

result  in  locally-reduced  flow  and  increased  deviation  angles  on  the 

suction  side  of  the  stator  blade  with  an  accumulation  of  more  slower 

moving  fluid  on  the  pressure  side.  The  deeper  stator  wake  regions  are 

those  corresponding  to  rotor  sampling  positions  YO  /Sr,  = 0.0,  0.69, 

and  0.83,  with  the  deepest  wake  region  being  associated  with  YO  /S  = 

K K 

0.83.  For  rotor  sampling  position,  YO  /S  = 0.83,  the  larger  tangential 
velocities  of  the  chopped  rotor  wake  fluid  appear  to  cancel  the 
smaller  tangential  velocities  of  the  stator  wake  fluid  with  the  result 
being  only  a small  amount  of  circumferential  variation  of  tangential 
velocities.  Also,  the  outward  radial  flow  of  the  chopped  rotor  wake 
fluid  cancels  the  inward  flow  of  the  stator  wake  fluid.  The  flow  angle 
tendencies  of  both  the  chopped  rotor  wake  and  the  stator  wake  seem  to 
add  together.  For  rotor  position  Y0„/Sn  = 0.5,  the  stator  wake  is 
least  affected  by  other  wake  fluid  as  It  is  most  like  an  Isolated 
stator  blade  wake.  Whenever  the  chopped  rotor  wake  is  clearly  discerni- 
ble at  the  measurement  plane,  Y0„/S,^  = 0.0,  0.69,  and  0.83,  it  is 
interacting  with  IGV  and  stator  wake  fluid,  and  thus,  its  distribution 
tends  to  be  characteristically  spread  out  and  shallow.  For  rotor 
position  YO^/Sj^  * 0.34,  little  of  the  chopped  rotor  wake  fluid  is  iden- 
tifiable at  the  measurement  plane. 

First  stator  exit  flow  data  for  10%  PHH  are  shown  in  Figure  13(b) 
with  related  cascade  plots  in  Figure  14(a).  As  the  rotor  blade  passes 
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In  stop-action  sequence  from  0.0  to  0.69,  the  effect  of  chopped 

rotor  and  IGV  wakes  on  the  stator  exit  flow  may  be  clearly  seen.  For 

rotor  sampling  position  YOj^/Sj^  • 0.0,  a nearly  isolated  (noninteracted) 

stator  wake  region  and  a separate  and  distinct  rotor  wake/lGV  wake 

interaction  region  are  present  at  the  measurement  plane.  For  Y0o/S„  • 

0.34,  the  stator  wake  region  is  slightly  deeper  with  a wake  sequence 

of  stator-rotor-IGV  for  Y/Sg  « 0.0  to  1.0.  The  rotor  wake/stator  wake 

interaction  occurs  just  upstream  of  the  measurement  plane,  and  thus, 

the  stator  wake  region  is  of  the  interacted  variety.  For  YO  /S  * 0.69, 

R R 

the  wake  sequence  is  stator-IGV-rotor , with  an  appreciably  different 

flow  pattern  from  the  ones  for  YO  /S  • 0.0  and  0.34.  For  Y0„/S„  * 0.69, 

R R R R 

the  large  axial  velocity  deficit  at  Y/Sg  = 0.6  is  the  result  of  an 

IGV  wake  region.  The  rotor  wake/stator  wake  interaction  occurs  downstream  j 

of  the  measurement  plane.  The  rotor  wakes  are  easily  identified  in  the  | 

tangential  angle  plots. 

Stator  exit  flow  data  for  30%  PHH  are  shown  in  Figure  13(c).  The 

deepest  stator  wake  region  occurs  for  rotor  sampling  position  Y0„/S  » 

R R 

0.69.  The  corresponding  cascade  plot  indicates  strong  stator  wake/ 

rotor  wake  interaction  at  the  measurement  plane.  The  least  influenced 

stator  wake  region  occurs  for  rotor  position  Y0_/S_  - 0.34.  In  this 

R K 

case,  the  velocity  and  angle  data  reflect  expected  trends  for  a 
noninteracted  stator  wake  region.  For  all  of  the  rotor  sampling 
positions,  the  tangential  angles  are  the  most  sensitive  Indicator  of 
rotor  wake  location  at  the  measurement  plane. 


Stator  exit  flow  data  for  70%  PHH  are  presented  in  Figure  13(d). 
The  rotor  wake  broadly  Influences  the  stator  exit  flow  and  makes  crisp 
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delineation  of  effects  difficult.  The  deepest  stator  wake  region  is 

measured  for  a rotor  sampling  position  of  YO  /S  = 0.0,  consistent 

R R 

with  the  cascade  plot  (see  Figure  14 [d]).  For  this  rotor  sampling 
position,  the  stator  wake/rotor  wake  interaction  is  occurring  just 
upstream  of  the  measurement  plane.  Radial  velocities  show  classical 
outward  rotor  wake  flow  and  inward  stator  and  IGV  wake  flow. 

Stator  exit  flow  data  for  90%  (Figure  13[e])  indicate  broad 
influence  of  the  rotor  on  the  stator  exit  flow.  Thus,  as  at  70%  span, 
wake  interaction  effects  are  relatively  difficult  to  sort  out.  Deeper 
stator  wake  regions  are  associated  with  stator  wake/rotor  wake  inter- 
action just  upstream  or  at  the  measurement  plane.  Tangential  angles 
are  again  good  indicators  of  chopped  rotor  wake  influence  location. 

By  looking  at  the  stator  exit  flow  from  hub  to  tip  for  three 
rotor  sampling  positions,  namely,  YO  /S„  = 0.0,  0.34,  and  0.69  (see 
Figure  13[f-h]),  several  observations  can  be  made.  Boundary  layer 
growth  has  extended  to  90%  PHH  behind  the  stator,  compared  to  only  95% 
PHH  behind  the  rotor.  An  upstream  effect  of  the  second  rotor  on  the 
flow  at  the  first  stator  exit  measurement  plane  is  not  noticeable, 
although  small  (less  than  4%)  upstream  effects  are  probably  present. 


5.  CONCLUSIONS 


Sunnnarized  below  are  the  conclusions  reached  after  studying  the 
periodic-average  flow  data  acquired  to  date.  Further  data  analysis 
and  acquisition  are  being  continued. 

The  IGV  suction  surface  flow  is  separated  at  all  radial  locations. 

The  separated  flow  appears  to  involve  unusual  unsteady  effects  that 
deserve  further  observation.  The  IGV  exit  flow  is  also  conplicated  by 
noticeable  first  rotor  upstream  effects  at  the  measurement  plane. 

The  chopped  IGV  wakes  move  downstream  within  fixed  avenues. 

Imbedded  rotor  and  stator  periodic  unsteadiness  was  found  to  be 
appreciable  in  portions  of  the  compressor  annulus,  depending  on  the 
extent  of  wake  interaction  involved  at  the  particular  location  considered. 
The  concept  of  chopped  wake  fluid  flowing  in  the  general  direction  of 
the  slip  velocity  appears  to  be  verified  by  the  present  data.  Chopped 
wake  interaction  with  the  suction  side  fluid  of  a blade  section  typi- 
cally resulted  in  locally  large  deviation  angles  and  small  axial 
velocities.  As  a result  of  wake  interaction,  tangential  and  radial 
velocities  were  observed  to  cancel,  while  tangential  angles  added. 

Behind  the  rotor,  large  differences  in  flow  fields  could  result  from 
rotor  wake/IGV  wake  interactions.  Behind  the  stator,  stator  wake/rotor 
wake  Interactions  did  not  lead  to  as  much  variation  in  flow  fields. 

Stator  wake/IGV  wake  interaction  effects  were  difficult  to  discern, 
even  though  IGV  effects  could  be  seen  at  most  spanwise  locations. 

Secondary  flow  effects  were  noticeable  in  some  cases  near  the  hub,  while 
leakage  and  scraping  effects  were  discernible  near  the  tip.  Outer 


78 


casing  boundary  layer  growth  was  evident  from  the  first  rotor  exit 
flow  to  the  first  stator  exit  flow. 

In  general,  the  periodic-average  flow  was  orderly  and  explainable 
with  consistent  wake  transport  and  interaction  plots  and  physical  | 

reasoning.  i 

! 
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I 7.  APPENDIX  A:  CALCULATOR  PROGRAMS  AND  STORAGE 


Calculator  data  acquisition  and  reduction  programs  used  in  the 
present  study  are  listed  in  this  section.  All  programs  and  data  are 
stored  on  cassette  tape  and  are  labeled  and  indexed  as  specified 
below: 

Flow  coefficient  program:  Calculation  of  overall  flow  coefficient  from 
flow  rate  venturi  meter  data.  Cassette  4C,  file  16. 

Actuator  position  correlation  program:  Linear  least  squares  correlation 
between  actuator  potentiometer  voltage  readout  and  actuator 
motion  for  probe  and  circumferential  positioning  actuators. 

Cassette  4C,  file  3. 


Hot-wire  effective  cooling  velocity /actual  velocity  ratio  calibration 
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8.  APPENDIX  B:  PARAMETER  EQUATIONS 

The  equations  used  in  the  periodic--average  measurement  system 
for  the  calibration  procedures  and  the  acquisition  and  reduction  of 
data  are  presented  below.  Symbols  and  notations  are  presented  on  page 
xii,  while  sign  conventions  are  generally  shown  in  Figure  10. 

8.1.  General  Parameters 


8.1.1.  Basic  Fluid  Properties 

2 

Barometric  pressure,  N/m  : 

P ^ = h,  [1.0-0.0018  (t,  -273.15)]  y, 

atm  hg@t,  baro  'hg@273  K 

baro 

3 

Density  of  air,  kg/m  : 

P 


(B-1) 


P = 


atm 
R t 


(B-2) 


Specific  weight  of  water,  N/m  : 


= ^ [ 996.86224  + 0.1768124  (j  t - 459.67) 


»2°  8, 

- 2.64966  X lO”^  4 " 459.67)^ 

+ 5.00063  X 10"^  " 459.67)^] 

8.1.2.  Blade-Element  Quantity 
Percent  passage  height  from  hub: 

8.1.3.  Miscellaneous 


(B-3) 


(B-4) 


Venturi  volume  flow  rate,  m /s; 


Q » 0.05229 

V 


^®c^H20^^vent 


(B-5) 
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Blade  velocity,  m/s: 


rirRPM 


(B-6) 


Calibration  nozzle  jet  velocity,  m/s: 


(B-7) 


Venturi  flow  coefficient: 


V A 


(B-8) 


8.2.  Three-Dimensional  Periodic-Average  Hot-Wire  Parameters 


Effective  cooling  velocity,  m/s: 

+ K2E|^  + K^E,^^  (B-9) 

Sensor  yaw  angle  relationship  (see  Figure  8): 

cos  a = cos  0-  cos  0 cos  0 + sin  0.  sin  0 (B-10) 

0 p y Op 

Effective  cooling  velocity /actual  velocity  ratio: 

V /V  = b.  + b,a  + b.0  + b,V  + b,a^  + b^0  ^ 

e 01  2p3  4 5p 

+ b,V^  + b_a6  + b-aV  + b.0  V (B-11) 

6 7 p 8 9 p 

Absolute  tangential  flow  angle  (see  Figures  8 and  9),  degrees: 

» B +0  + 0 (B-12) 

0 rav  a, off  y,c 

Radial  flow  angle  (see  Figure  8),  degrees: 


6 * -0 
r p 


(B-13) 


Radial  component  of  fluid  velocity,  m/s: 


V - V sin  B 
r r 


(B-14) 
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Axial  component  of  fluid  velocity,  m/s: 

V = V cos  6 cos  6„ 
z r 0 

Tangential  component  of  absolute  fluid  velocity,  m/s: 
Vq  = V cos  6^  sin  0^ 

Tangential  component  of  relative  fluid  velocity,  m/s: 


Relative  fluid  velocity,  m/s: 

V'  + (V^)^ 

Relative  tangential  flow  angle,  degrees: 
6*  = sin"l  (V-  /V,) 


(B-15) 


(B-16) 


(B-17) 


(B-18) 


(B-19) 


9.  APPENDIX  C:  TABULATION  OF  PERIODIC-AVERAGE  HOT-WIRE  DATA 


The  periodic-average  circumferential  survey  data  are  tabulated  in 
this  section.  The  data  represent  flow  field  parameters  downstream  from 
the  IGV  row  (station  2),  first  rotor  row  (station  3),  and  the  first 
stator  row  (station  4).  Each  vector  is  completely  defined  by  velocity 
magnitude,  tangential  flow  angle,  and  radial  flow  angle.  The  computer 
headings  are  defined  as  follows: 

Y/SS  = circumferential  spacing,  Y/S 
V = absolute  velocity,  V,  m/s 

BETA  Y = absolute  tangential  flow  angle,  gO,  degrees 
BETA  R = radial  flow  angle,  g^,  degrees 
PHH  * percent  passage  height  from  hub,  PHH 

YOR/SR  = circumferential  ratio  blade  sampling  position,  YOn/S_ 

K R 


Table  C-1.  Hot  wire  circumferential  survey  data  obtained  with  the  periodic-average  measure 
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